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We present a search for B decays to a charged scalar meson and a ir° where the decays 
to an r\ meson and a n + . The analysis was performed on a data sample consisting of 383x10 BB 
pairs collected with the BABAR detector at the PEP-II asymmetric-energy B Factory at SLAC. We 
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find no significant signal and set an upper limit on the product branching fraction 
B(B + -> a^iv ) x B(a£ -> r?7r + ) of 1.4xl0" 6 at the 90% confidence level. 

PACS numbers: 13.25.Hw, 12.39.Mk 



The structure of scalar mesons is a subject of some 
debate 0, Q. Proposed models include two-quark or 
four-quark states with potential contributions from glue- 
balls or a molecular admixture of KK meson pairs. 
Measurement of the branching fraction for the mode 
B + — > ciqIt Q is expected to provide an effective test 
of the two- and four-quark models [f| . The Feynman di- 
agrams for the decay in the two-quark case are shown in 
Figure [T] Those for the four-quark case are similar ex- 
cept for an ss pair produced from the vacuum internal 
to the <Iq meson. The color-allowed electroweak tree dia- 
gram shown in Figure^a) is suppressed for all models 
since the W + is constrained to decay to states of even G- 
parity (a generalization of C symmetry to cover particle 
multiplets) within the Standard Model, whereas the 
has odd G-parity [1| . This diagram is also suppressed due 
to vector current conservation considerations. Therefore, 
the color-suppressed tree diagram in Figure 0Jb) and the 
helicity-suppressed electroweak annihilation diagram in 
Figure [TJc) become important. The gluonic penguin pro- 
cess in Figure [ljd) is highly suppressed and is therefore 
not expected to contribute significantly. 

The amplitudes for the above diagrams depend on the 
model used; in particular the annihilation diagram 
is heavily suppressed in a four-quark model. Hence mea- 
surement of the branching fraction provides the potential 
for model discrimination. In the two-quark case, thepre- 
dicted branching fractions go as high as 2xl0~ 7 [5j, |f|. 
However, in the four-quark case the prediction for the 
branching fraction is an order of magnitude lower. 

The branching fraction for the result quoted below 
will be given in terms of the product B{B + — > a^n ) x 
B(cIq — > rjir + ) since the branching fraction B(a^ — > 7jtt + ) 
is not well measured, although it is thought to be ap- 
proximately 85% [l[. 

The analysis presented in this paper is based on 
347 ft)" 1 of data collected at the T(4S) resonance with 
the BaBAR detector at the PEP-II asymmetric-energy 
e + e~ collider located at the Stanford Linear Accelera- 
tor Center. This corresponds to (383±4)xl0 6 BB pairs. 

The BABAR detector has been described in detail previ- 
ously [~7j . Track parameters of charged particles are mea- 
sured by a combination of a 5-layer double-sided silicon 
vertex tracker and a 40-layer drift chamber (DCH), both 
operating in the 1.5 T magnetic field of a superconduct- 
ing solenoid. Photons and electrons are identified using 
a CsI(Tl) electromagnetic calorimeter. Further charged 
particle identification (PID) is provided by measurements 
of the average energy loss (dE/dx) in the tracking devices 
and by an internally-reflecting, ring-imaging Cercnkov 
detector (DIRC) covering the central region. 




u ► u u ► u 

(a) (b) 




(c) (d) 



FIG. 1: The Feynman diagrams contributing to the process 
B + — > (IqTv in the two-quark model, (a) is the external 
(color-allowed) tree, (b) the internal (color-suppressed) tree, 
(c) the annihilation process and (d) the gluonic penguin pro- 
cess. 



The analysis focuses on Oq mesons produced from the 
decay B + — > ciqTt , followed by a q — > ?77r + , where the r\ 
meson subsequently decays to 77 or 7r + 7r~7r° final states. 
The 7T° mesons used are reconstructed via the decay 
7T° — > 77. The selections used for the analysis are the 
result of an optimization procedure based on ensemble 
Monte Carlo (MC) studies. In these studies, a sample of 
MC candidates is produced for given selection criteria by 
generating randomly from probability density function 
(PDF) distributions defined with the selection applied. 
By re-fitting to the datasets for each set of selection cri- 
teria it is possible to select the set that yields the max- 
imum sensitivity to signal. This is done independently 
for each decay mode considered. In both cases candi- 
dates are required to satisfy 0.8 < < 1.2 GeV/c 2 with 
the rj candidates satisfying 0.51 < m 77 < 0.57GeV/c 2 or 
0.540 < m 37r < 0.555 GeV/c 2 . The ir° produced from 
the 77 — ► 7r + 7T~7r° decay is required to satisfy 0.10 < 
m^o < 0.16GeV/c 2 . The ir° daughter of the B candidate 
is required to satisfy 0.115 < m^o < 0.150 GeV/c 2 . This 
selection is tighter than for the 7r° produced from the 77 
meson since it is of significantly higher energy and there- 
fore has a better resolution. The charged track from the 
candidate decay is required not to satisfy DIRC- and 
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DCH-based PID criteria consistent with a kaon hypoth- 
esis. This PID selection has been measured to be more 
than 80% efficient for tracks with momenta up to 4 GeV/c 
with a pion mis-identification rate lower than 10% over 
the same range. 

A B meson candidate is characterized kinematically 
by the energy difference AE = Eb — \\fs and energy- 
substituted mass mEs = {\s — pfj) 5 > where s is the 
square of the centre-of-mass energy of the colliding 
beams, (Eb,Pb) is the candidate B meson 4-momentum 
and all values are expressed in the T(4S) frame. Signal 
events peak around zero for AE, and at the B meson 
mass for rngs- The resolutions for AE and Toes are ap- 
proximately 30MeV and 3MeV/c 2 , respectively. We re- 
quire \AE\ < 0.35 GeV and 5.20 < to es < 5.29 GeV/c 2 as 
an input for the fit used to extract signal and background 
parameters (described below) in order to maximize the 
available statistics. 

The principal source of background in the analysis 
arises from random combinations in continuum e + e~ — > 
QQ (Q = M , d, s, c) events. These contributions are reduced 
in part by placing a selection on the variable | cos(#tb)|, 
where Qtb is the angle between the thrust axis of the 
B candidate and the thrust axis of the rest of the event 
calculated in the T(AS) frame. Candidates formed in jet- 
like qq events will peak at | cos(#tb)| values approaching 
1, whereas signal B decays will follow an almost flat dis- 
tribution as they are isotropic in this angle. We require 
|cos(6>tb)| < 0.7 for both r\ channels. The final vari- 
able used in the analysis is a linear Fisher discriminant 
T that consists of the angles of the B momentum and 
B thrust axis (in the T(AS) frame) with respect to the 
beam axis, and the zeroth and second Legendre moments 
of the energy flow computed with respect to the B thrust 
axis § . The reconstruction efficiencies after selection are 
presented in Table Q] 

The analysis uses an extended unbinned maximum- 
likelihood fit to extract yields for the modes under study. 
The input variables to the fit are AE, toes, F an d the 



candidate resonance mass 
function for the fit is defined as: 



The extended likelihood 



C = 



N\ 



1 N 

-n 

i=l 



M 



(1) 



where Vj is the normalized PDF for a given fit component 
j. For each candidate % the PDF is evaluated using the 
fit variables of that candidate. The M fit components 
are the signal and all background contributions. The 
total number of candidates is given by N with the yield 
associated with each fit component given by rij . The fit 
for each r\ channel consists of 16 components modeling 
signal and continuum candidates separately as well as 
charged and neutral charmed B meson decays. There 
are then 12 components modeling individual charmless 
modes which were found to contribute a background to 



the signal. The yields for all B background components 
are held fixed in the final fit using values calculated from 
the latest branching fraction estimates whereas the 
signal and continuum background yields are allowed to 
vary. 

The fit model is constructed in order to extract signal 
candidates effectively from a sample where multiple re- 
construction hypotheses exist for each event. The signal 
MC events have an average candidate multiplicity of 1.4 
for both r\ decay modes. 

In this analysis separate PDFs were used to discrimi- 
nate between correctly and incorrectly reconstructed sig- 
nal candidates in MC. This was achieved by using MC 
information to separate the signal MC candidates into 
an almost pure sample of correctly reconstructed candi- 
dates and a sample consisting mainly of incorrectly re- 
constructed candidates. By iteratively fitting the sepa- 
rate PDFs to each sample in turn, a consistent set of 
PDFs for the two cases was obtained. The component 
for correctly reconstructed candidates was then taken to 
model signal candidates in the final fit to data. The frac- 
tion of events in the MC that were identified as correctly 
reconstructed by the fit was approximately 62% for both 
r\ channels. The signal candidate yield resulting from the 
fit to MC was verified to be consistent with that expected. 

The shapes of the distributions for incorrectly recon- 
structed signal were found to be similar to continuum 
background and thus any such candidates are assumed 
to be absorbed into the yield associated with the contin- 
uum PDF. Modeling signal candidates in this way was 
shown using ensemble MC studies to provide better sen- 
sitivity to signal than other methods. As a final test, 
the method was validated using ensemble MC studies to 
show that it introduced no bias into the final fit result. 

Any continuum and BB backgrounds that remain after 
the event selection criteria have been applied are identi- 
fied and modeled using Monte Carlo simulation based on 
the full physics and detector models [lOj. Charmless B 
decays providing a background to the signal are identi- 
fied by analyzing the MC candidates passing selection 
from a large mixed sample of Standard Model B decays. 
Charged and neutral charmed B decays are modeled sep- 
arately and individual components are included for each 
charmless B decay mode found to contribute. The PDF 
parameters for each B background component are ob- 
tained from MC samples and held fixed in the final fit to 
data. Those for the continuum background shape are left 
free in the final fit. The contributions from two charmless 
backgrounds with the same final state as signal, those for 
B + —* a (1450) + 7r° and non-resonant B + — > rjTr + ir°, are 
estimated using fits to the relevant regions of the Dalitz 
plane. Any potential interference effects were neglected 
since the fits gave no significant yields for these modes. 

The total PDFs are modeled as products of the PDFs 
for each of the four fit variables. The signal shapes 
in AE, toes, tn^ix and T are modeled with a Novosi- 
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TABLE I: The results of the fit to the full data set, and other values required for calculating the branching fraction. All B 
background yields were held fixed. The upper limit is shown first with only the statistical error and then with the total error. 



Required Quantity/Result 


j) -> 77 


Tj — » 7T + 7T 7T U 


Candidates to fit 

Fixed B Background (candidates) 
Signal Yield (candidates) 
Continuum Yield (candidates) 
ML Fit Bias (candidates) 


f03054 

1640 
-8 ± 19 
101400±300 
5.2±3.0 


31626 

942 
13±13 
30700±200 
-2.0±1.3 


Efficiencies and BFs 






Efficiency (%) 
B(r, -> X) (%) 
Branching Fraction (xlO 6 ) 


16.3±0.1 
39.4±0.3 
-0.6«; 8 7 (stat) t° .t (syst) 


10.2±0.1 
22.6±0.4 
\.7+_\l (stat) t° .l (syst) 


Combined Mode Results 




Branching Fraction (xl0~ 6 ) 
Significance 

Upper Limit 90% C.L. (xf(T 6 ) 
Upper Limit 90% C.L. (xlO -6 ) 


0.11°;? (stat) ±° 3 3 (syst) 

0.1a (stat + syst) 
< 1.3 (statistical error only) 
< 1.4 (total error) 



birsk function, the sum of two independent Gaus- 
sians, a Breit-Wigner, and an asymmetric Gaussian, re- 
spectively. The signal parameters used for the line- 
shape are a Breit-Wigner peak value of 983 MeV/c 2 with 
a width of 79MeV/c 2 . These were used in the MC sim- 
ulation and are consistent with previous analyses [T3 |. 
although the width is considered to be uncertain over 
a conservative range of 50-100 MeV/c 2 in the evaluation 
of systematic error. Slowly-varying background distribu- 
tions in T and are modeled with Chebychev poly- 
nomials of the appropriate order. Such polynomials are 
also used for AE in the charmed B and continuum back- 
ground cases. For these components tties is modeled 
with an ARGUS [l3| threshold function. In the case of 
charmless B backgrounds, AE and 771.es are modeled 2- 
dimensionally using non-parametric PDFs so as to 
model correlations between the two variables. Studies of 
the MC samples for each mode have shown that these 
correlations can be as high as 29%. 

The results of the analysis are presented in Table Q] 
The statistical errors on the signal yields are defined us- 
ing the change in the central value when the quantity 
— 21n£ increases by one unit from the minimum. The 
significance is taken as the square root of the difference 
between the value of —2 In £ for zero signal and the value 
at the minimum (including additive systematics) . 

For the purposes of the branching fraction calculation 
we assume that the T(AS) decays with an equal rate to 
both B + B~ and B°B° 15|. The fit bias is measured us- 
ing an ensemble MC study based on a parameterization 
taken from the fit to data with all yield values taken from 
data. Where a negative yield is found a value of zero is 
used for the study. The branching fraction results from 
the two 77 decay modes are combined by forming the prod- 
uct of the likelihood functions, after their maxima have 



been shifted to account for fit bias. The functions them- 
selves are defined by computing the likelihood values for 
signal yields around the maximum. Systematic errors are 
included at the required stages in the calculation depend- 
ing on correlations between the two rj channels. 

We find no significant signal in either r] decay mode and 
thus quote upper limits on the branching fraction at the 
90% confidence level (C.L.), taken to be the branching 
fraction below which lies 90% of the total of the likelihood 
integral in the positive branching fraction region. 

In Figure [5] we show projections of each of the four fit 
variables for both the r\ — > 77 and 77 — * 7r + 7r~7r° decay 
modes. To enhance the visibility of a potential signal, the 
candidates in these figures have been required to satisfy 
the condition that the likelihood ratio C S i g /[C S i g +S Cbkg] 
for any candidate be greater than 0.6. Here Cx is the 
likelihood for a given event being described by either the 
signal or background model. The likelihoods are calcu- 
lated for each figure separately, excluding the variable 
being plotted. As can be seen there is no significant sig- 
nal peak for either mode. 

The largest sources of systematic uncertainty in the 
analysis arise from poor knowledge of the lineshape 
and from the error in the estimated background contri- 
butions. By varying the width of the Breit-Wigner 
between 50 and 100MeV/c 2 we predict an uncertainty of 
approximately +5 and —4 candidates for rj — > 77 and 
+0.5 and —1 candidate for 77 — > n + ir~iT Q . Varying the 
charmless yields within their branching fraction errors (or 
±100% where a limit is used), and the charmed B yields 
by ±10%, gives an estimated uncertainty of ±4 candi- 
dates in r] —> 77 and ±1 candidate in 7; — > 7r + 7r~7r . The 
error due to the uncertainty in the fit bias was calculated 
as the sum in quadrature of 50% of the measured bias 
and its statistical error, as taken from the ensemble MC 
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FIG. 2: Likelihood-ratio-enhanced projections for the four fit variables (left to right) for the r\ — > 77 (top) and 77 — > 7t + 7t~7t 
(bottom) cases. Experimental data are represented by points with error bars, solid blue curves represent the full fit model. For 
the -qz-n case, the combined background component is represented by the black dash-dotted curve and the signal component by 
the red dashed curve. The efficiency of the likelihood ratio selection on the signal component in the 773^ case is 80.8%. 



TABLE II: Estimated systematic errors in the final fit result. 
Error sources which are Correlated and Uncorrelated for the 
two 77 decay modes are denoted by [C] and [U], respectively. 



Source of Uncertainty 


77 -> 77 


77 —> 7T + 7T 7T U 


Additive (Candidates) 






Fit Parameters [U] 


+5.9 
-4.4 


+ 0.5 
-1.8 


Charmless Yields [U] 


+3.6 
-3.7 


+ 1.2 
-1.2 


Charm Yields [U] 


+0.2 
-0.3 


+0.2 
-0.2 


Fit Bias [U] 


±3.0 


±1.3 


Total Additive (Candidates) 


+7.5 
-6.5 


+1.9 
-2.6 


Multiplicative (%) 






Neutral Efficiency [C] 


±6.0 


±6.0 


Tracking Efficiency [C] 


±0.5 


±1.4 


|cos(#tb)| Selection [C] 


±3.0 


±3.0 


MC Statistics [U] 


±0.4 


±0.3 


Number of BB Events [C] 


±1.1 


±1.1 


Daughter r\ Decay BF [U] 


±0.7 


±1.8 


Total Multiplicative (%) 


±6.9 


±7.2 


Total BF Syst Error (xlCT 6 ) 


+0.4 
-0.3 


+0.3 
-0.4 



study described above. This value was calculated to be 
approximately ±3 candidates in the 77 — > 77 channel and 
±1 candidate for 77 — > 7r + 7r~7r°. 

Further sources of systematic uncertainty, which are 
multiplicative rather than additive, affect the efficiency 
and thus enter into the branching fraction calculation. 
Limited signal MC statistics account for 0.4% in both 77 
decaymodes. Auxiliary studies on inclusive control sam- 
ples [8], predict errors of 0.5% per charged track and 3%. 
per reconstructed 77 or 7r° decaying to two photons. The 
estimate of the number of produced BB events is uncer- 



tain by 1.1%. The uncertainties in B daughter product 
branching fractions are taken to be 2% for 77 — > 77 and 
3% for 77 — > 7r + 7r~7T° A summary of all systematic 
error contributions is presented in Table |TTJ 

In conclusion, we do not find a significant signal for the 
mode B + — > a^ir . We set an upper limit at 90% C.L. on 
the branching fraction B(B + — > a^zr ) x B(a^ — > 777r + ) 
of 1.4xl0 -6 , suggesting that there is insufficient sensi- 
tivity with the current dataset to probe the predicted 
theoretical parameter space, with the largest predicted 
branching fraction being 2xl0~ 7 [5]. We are therefore 
unable to comment on the validity of any of the current 
models of the . 

We are grateful for the excellent luminosity and ma- 
chine conditions provided by our PEP-II colleagues, and 
for the substantial dedicated effort from the comput- 
ing organizations that support BABAR. The collaborat- 
ing institutions wish to thank SLAC for its support and 
kind hospitality. This work is supported by DOE and 
NSF (USA), NSERC (Canada), IHEP (China), CEA and 
CNRS-IN2P3 (France), BMBF and DFG (Germany), 
INFN (Italy), FOM (The Netherlands), NFR (Norway), 
MIST (Russia), MEC (Spain), and STFC (United King- 
dom). Individuals have received support from the Marie 
Curie EIF (European Union) and the A. P. Sloan Foun- 
dation. 



* Deceased 

f Now at Tel Aviv University, Tel Aviv, 69978, Israel 

* Also with Universita di Perugia, Dipartimento di Fisica, 



8 



Perugia, Italy 
® Also with Universita della Basilicata, Potenza, Italy 
■ Also with Universitat de Barcelona, Facultat de Fisica, 
Departament ECM, E-08028 Barcelona, Spain 
[1] Particle Data Group, Y.-M. Yao et al, J. Phys. G33, 1 

(2006), "Note on Scalar Mesons", p. 546. 
[2] V. Baru et al, Phys. Lett. B 586, 53 (2004). 
[3] Throughout this paper, charged-conjugate decays are 
also implied. 

[4] S. Laplace and V. Shelkov, Eur. Phys. Jour. C 22, 431 
(2001). 

[5] D. Delepine, et al, Eur. Phys. Jour. C 45, 693 (2006). 
[6] H.-Y. Cheng, et al, Phys. Rev. D 73, 014017 (2006). 
[7] BABAR Collaboration, B. Aubert et al, Nucl. Instrum. 

Methods Phys. Res., Sect. A 479, 1 (2002). 
[8] BABAR Collaboration, B. Aubert et al., Phys. Rev. D 70, 

032006 (2004). 

[9] Particle Data Group, Y.-M. Yao et al., J. Phys. G33, 
1 (2006); Heavy Flavour Averaging Group (HFAG), E. 
Barberio, et al, hep-ex/0603003 (2006). 



[10] The BABAR detector Monte Carlo simulation is based on 
GEANT4: S. Agostinelli et al, Nucl. Instrum. Methods 
Phys. Res., Sect. A 506, 250 (2003). 

[11] The Novosibirsk function is defined as 

f(x) = A s exp(-0.5(ln 2 [l + At(x - x )]/t 2 + r)) where 
A — sinh(r Vln 4)/(<tt Vhi 4) , the peak is xo, r is the tail 
parameter and A 3 is a normalization factor. 

[12] E852 Collaboration, S. Teige et al., Phys. Rev. D 59, 
012001 (1998). 

[13] ARGUS Collaboration, H. Albrecht, et al, Phys. Lett. 
B 241, 278 (1990); Function defined as 
f(x) = Nx[(l-(x/E beam )) 2 exp\p(l-(x/E bearn )) 2 ]] 1/2 
where iV is a normalization factor, p a shape parameter 
and Eteam is 50% of the centre-of-mass energy of the 
colliding beams. 

[14] K. S. Cranmer, Comput. Phys. Commun. 136, 198 
(2001). 

[15] See for instance BABAR Collaboration, B. Aubert et al., 
Phys. Rev. D 69, 071101 (2004), and references therein. 



